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ABSTRACT 

The balloon-borne Solar Disk Sextant (SDS) experiment has measured the 
angular size of the Sun on seven occasions spanning the years 1992 to 2011. The 
solar half-diameter - observed in a 100- nm wide passband centered at 615 nm - 
is found to vary over that period by up to 200 mas, while the typical estimated 
uncertainty of each measure is 20 mas. The diameter variation appears as if it 
might be cyclic, although it is not in phase with the solar activity cycle; thus, 
the measured diameter variation is not simply an observational artifact of surface 
activity. The SDS is described here in detail, as is the complete analysis procedure 
necessary to calibrate the instrument and allow comparison of diameter measures 
across decades. 

Subject headings: Sun: fundamental parameters — astrometry — balloons 

1. Introduction 

Measuring the solar diameter (and its variations) has a long and controversial history 
with results highly inconsistent with each other, not only in absolute values, but also in 
trends with time. Detailed reviews of these efforts are presented in Djafer et al. (2008) and 
Thuillier et al. (2011) and will not be repeated here. More specifically, here we will only 
address direct measurements, as opposed to values derived from the analysis of historical data 
(e.g. duration of solar eclipses, timings of transits of Mercury, etc.). The main conclusion 
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of the paper by Djafer et al. was that measuring the solar diameter is a very difficult 
undertaking, and that the principal causes of the inconsistent results are the effects of the 
Earth's atmosphere, differences in the definition of the limb edge, the spectral range of the 
measurements, and instrumental variations that cannot be independently calibrated with a 
precision higher than the expected variations. 

Whereas differences of the results between the different diameter determinations can 
arise from the various edge definition algorithms, and the range of wavelength of measure- 
ments, the extreme differences in trends can only arise from atmospheric effects, and/or 
from the large instrumental effects produced by the extreme environment of Sun-pointing 
telescopes that do not have a system for internal scale calibration. 

If we confine our discussion to space borne (or space-like) measurements, the atmo- 
spheric effects are removed, or greatly reduced, and what remains are the instrumental ef- 
fects. This selection leaves only 5 experiments for which diameter measurements have been 
attempted: the MDI on SOHO, the HMI on SDO, RHESSI, the SODISM experiment on 
PICARD, and the SDS. Of those, only the SDS and SODISM have internal calibration capa- 
bility. Because of the absence of such calibration, the RHESSI results only presented determi- 
nations of the solar oblateness, but not of the diameter (Fivian et al. 2008). This was based 
on the reasonable assumption that the instrumental scale would not change significantly 
during the short interval in which an instrumental rotation is performed. Angular diameter 
determinations with SOHO/MDI (Bush et al. 2010) find no discernible variation, although 
without an internal angular calibration these measurements must rely on other means of 
correcting significant long-term changes in the instrument. As with RHESSI, SDO/HMI has 
made observations sufficient to measure the solar oblateness (Fivian et al. 2012), although 
results of these measures have yet to appear. SDO has not been operating long enough 
to address the issue of long-term solar diameter variation. PICARD/SODISM (Thuillier et 
al. 2008, Assus et al. 2008) is still in a validation phase, and no results have been published to 
date. As a consequence, the only currently available results from instruments in a space-like 
environment and with internal calibration, are the SDS results presented here. 

Because the SDS is balloon-borne, it can only be flown during periods when the high- 
velocity stratospheric winds change direction (twice per year) at which time, and for a few 
days, their velocity is low (the so called turnaround period). During turnaround we can 
have day-long flights and still remain within the range of control of the NASA/Columbia 
Scientific Balloon Facility (CSBF), the organization that operates these flights. Moreover, 
at the location of the CSBF base in Fort Sumner, New Mexico, turnaround occurs in early 
May and late September. At the spring turnaround, the Sun is too high in the sky, so 
that for several hours around local noon the Sun would be behind the balloon, and hours 
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of observation would be lost. As a consequence, flights are not recommended at that time. 
Hence, only fall flights are desirable. More precisely, since the SDS in the current gondola 
cannot tilt above 53°, those flights have to take place on September 25 or later. 

Although we have flown the SDS eleven times, starting in the late 1980s, the instru- 
mental configuration required to reach a precision of tens of milliarc sec (to be explained in 
Section [2]) was only achieved starting with Flight 6, in 1992. As a consequence, this paper 
only presents the results of the seven flights since and including Flight 6. Of those, Flight 10 
has a significantly small number of observations because of a failure of the recording device 
shortly after the balloon reached float altitude. The flights with the appropriate precision 
were carried out on the dates shown in Table [2j 

This paper describes the SDS instrument and its operation in Sections [2] and [31 re- 
spectively. Section H] details the data-analysis procedures and pipeline. Section presents 
the diameter results and Section [6] discusses them, while Section [7] provides a concluding 
summary. 

2. SDS Instrument 

The current version of the SDS balloon borne instrument is a package that has been 
successfully flown several times on 12 - 29 Mcf (million cubic ft) balloons. The basic principle 
of the SDS instrument is the use of a mechanically and optically stable beam splitting wedge 
(BSW) as an angle reference to form a double image of the Sun separated by slightly more 
than its angular diameter. The constancy of the angle of the BSW is achieved by utilizing 
molecular contact fabrication techniques. It can be shown that by measuring the small 
distance between images, d, one can achieve the necessary accuracy much more easily than 
if one attempts to measure the full diameter directly. The level of dimensional stability 
required within the focal plane is relaxed as d/D, where D is the distance between the 
center of two solar images; in our case d ~ 6 mm and D ~ 200 mm. Figure [TJ illustrates the 
measurement concept and the layout of the detectors in the SDS focal plane. 

This method also allows one to monitor and correct for instrument changes (e.g. focal 
length changes, etc.) in a straightforward way. Changes in any of the optical components 
downstream of the wedge will affect both images and can be calibrated out. The instrument 
consists of three basic items, the optical system, the assembly of linear-array detectors, and 
the electronics package mounted separately from the telescope. The accuracy of the SDS 
derives from the system design, which uses a single optical train to transfer the split solar 
images to the detectors. 
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Fig. 1. — Upper: Optical layout of the SDS. The SDS is built around a ruggedized Questar 
telescope supplemented with a beam splitting wedge and a Barlow relay lens to increase 
the focal length. The beam splitting wedge produces a direct image and a reflected image 
that is offset by slightly more than the Sun's diameter. The solar limb is measured along ten 
radials by seven linear-array CCDs mounted in the focal plane. Insert: Detail of the detector 
plane. The red line segments indicate the location of the CCDs; the dotted circle shows the 
approximate field of view. The solar diameter is determined by precise measurement of the 
disks' separation, D, and gap, d. 
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In the current instrument, F = 20.5 m, and 2W = 1978.94". Further details of the 
instrumental properties are given in Sofia et al. (1991) and summarized in Table [H 

2.1. Optical system 

The stability of the optics needed by this instrument is of the level commonly used 
in optical interferometry. For this reason, the balloon flight version currently utilizes simi- 
lar techniques (molecular bonding) and materials (quartz and Zerodur). Thus, the optical 
system consists of the following components: 

1. beam splitting wedge (BSW), 

2. Cassegrain telescope, 

3. relay lens, to achieve the required focal length, and 

4. detector support. 

2.1.1. Beam splitting wedge (BSW) 

Because the wedge is the most critical element in the optical system, great care is taken 
with its design and manufacture. The wedge consists of two fused silica flats separated by 
an annular silica ring polished to an angle of about 1000 arcsec. Molecular contact bonding 
is used to hold the assembly in alignment. The surfaces are flat to 1/50 wave at 630 nm 
and have dielectric coatings to define the bandpass and reduce the solar transmission to an 
acceptable level. The mirrored surfaces have a high reflectivity (> 0.9) so that the intensities 
of the two images are approximately in the ratio of 5:4. The passband transmitted by the 
series of coatings on the wedge surfaces is centered at approximately 615 nm and is roughly 
100 nm wide. 



Table 1: SDS characteristics 



focal length : 


20.5 m 


aperture : 


12.7 cm 


passband : 


A = 615 nm, SX ~ 100 nm 


pixel scale : 


0.13 arcsec/pixel 


wedge angle x2 : 


1978.94 arcsec 
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2.1.2. Telescope 

The 17.8 cm aperture ruggedized Questar telescope is used with a reduced aperture of 
12.7 cm to accommodate the effective wedge aperture, and has a nominal focal length of 2.5 
m. This instrument has optical elements made of fused silica and Zerodur, and a main body 
and mounting parts made of Invar to minimize thermal effects. 

2.1.3. Relay lens 

In order to provide the appropriate plate scale to match the resolution of available 
Charge Coupled Device (CCD) detector elements, a magnification of the image is required. 
A magnification by a factor of 8, to an effective focal length of 20 m, is provided by a 
multi-element Barlow lens. 



2.1.4. Detectors 

The linear-array CCDs are mounted on a ceramic printed circuit board having a low 
thermal coefficient of expansion near that of the CCD cases and a high thermal conductivity. 
This provides uniformity of expansion for the entire assembly during temperature excursions, 
with the amplitude of the excursions being both small and well characterized. The positions 
of the endpoints of the CCD elements are determined by measurement on the Yale University 
microdensitometer to a precision of 2 microns. This assembly of detectors is held in the focal 
plane of the instrument. Detection of the solar images is achieved by use of seven such CCD 
arrays. The central detectors are used to measure the gap, while the outer detectors are 
used, in conjunction with those in the center, to define the centers of the images and thus 
the plate scale. The linear arrays are Texas Instrument virtual phase 1728 element devices 
having pixel dimensions of 12.7 x 12.7 microns and packaged in a narrow, windowed carrier. 
The CCD array is calibrated in the laboratory using an integrating sphere as a uniform light 
source. A bandpass filter provides the proper wavelength of light to match the solar input. A 
geometry representative of that in the flight optical system with regard to source diameter, 
angle subtended, etc. is used. A zero-level offset and linear and quadratic sensitivity factors 
are determined for each of the pixels in the seven CCDs. This calibration procedure has 
produced an effective pixel uniformity of better than 1/4%, which is easily adequate for our 
present data analysis procedures. In order to further improve the calibration, during flight 
we stop Sun pointing at a given time, and let the pointing drift over the solar disk. 
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2.2. Pointing System 

The analog pointing system, which utilizes both a LISS (Lockheed Intermediate Sun 
Sensor) and real-time feedback from the detectors, achieves an operational stability equal to 
or less than 12 arcsec over the entire flight (10-12 hours). 

2.3. Analog data system 

The function of the SDS analog data system is to control data transfer from the CCDs to 
buffers that can be read by the digital data system. Sample/Hold buffers are used to freeze 
the pixel signal and to set the base level equal to the black reference level for 12 bit analog 
to digital conversion. Each of the CCDs has its own system and they are read in parallel to 
freeze the image. These conversions are performed at the (commandable) exposure rate with 
the digitized values serially shifted to the 7 memory buffers in the onboard computer at up 
to 4 Mbits/sec/channel. The analog data system is mounted near the detectors to minimize 
noise. 



2.4. Digital data system 

The digital data system is based on a 66 MHz 32-bit Intel 486 processor together with 4 
Mbytes of main memory. It uses a pair of 4 Gbyte disk drives for data and program storage, 
and 7 - 64Kbyte data RAM buffers. Seven channels of serial data at up to 4Mbits/sec/channel 
are used to fill the 7 - 64Kbyte buffers. A CCD clock pulse is generated at 4 times the pixel 
rate which is command programmable for a rate of 4 to 20 microsec/pixel. Detector power 
and control signals are supplied by the computer. Functions of the computer are as follows: 

1. provide a command link between the SIP (Standard Instrument Package) and the SDS 
instrument, 

2. provide a data link to the ground via the SIP, 

3. provide a commandable CCD pixel clock, 

4. receive 7 channels of serial pixel data along with a bit and frame sync, 

5. store 18 digitized data frames for each of 7 CCDs, (which we refer to as a data cycle), 

6. calculate edge data and supply computed pointing position corrections to the onboard 
computer for fine pointing, 

7. when pointing accuracy is satisfied, calibrate the data sets and compute solar edge 
positions, 



-8 - 



8. periodically obtain a set of data from the 7 CCDs for transmission to the ground. 

Steps vi) and vii) refer to in-flight determinations of the solar edges for the purpose of 
fine tuning the pointing. Section H] describes the much more thorough edge-determination 
analysis performed on the data after the flight. 

3. Observations 

Table [2] lists the dates of the seven SDS balloon flights that have yielded useful mea- 
surements of the solar diameter. All flights were conducted by the CSBF from their station 
in Fort Sumner, New Mexico. The table lists peak altitude for each flight, which is near the 
float altitude at which the observations are made. 

Flight 6 was the first to be made with the contact-bonded wedge. A spring mechanism 
had been used in previous flights to try to maintain the geometry of the wedge but this proved 
to be totally inadequate, thus rendering data from the earlier flights unusable. Between 
Flights 6 and 7, the SDS was sent to the White Sands Missile Range optical instrumentation 
laboratory for a thorough refurbishment; this included making structural changes to the 
primary mirror mount and locking the previously adjustable focus in place following a careful 
benchtop alignment. A failure of the onboard hard drive during Flight 10 meant that only 
a small fraction of the in-flight data was recorded, as part of the telemetry stream sent to 
the ground station. Still, there is a sufficient amount of data from Flight 10 to allow for a 
worthwhile radius measurement to be made. 

Flight 11 employed a substantially larger balloon than the other flights, (29 Mcf versus 
12 Mcf), accounting for its somewhat higher peak altitude. An examination of the width of 



Table 2: SDS flights 



Fit. 


Date 


Alt. 


Comments 


# 




(km) 




6 


1992-09-30 


30.7 


First flight using contact-bonded wedge 


7 


1994-09-26 


31.6 


First flight with fixed focus 


8 


1995-10-01 


31.5 




9 


1996-10-10 


32.0 




10 


2001-10-04 


31.6 


Onboard memory failure; < 1% of data retrieved 


11 


2009-10-17 


33.3 


Anomalous excursion in image quality 


12 


2011-10-15 


31.0 
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the solar limb edge (as determined using the procedure described in Sect H|) over the course 
of the flight indicates an anomalous behavior. Figure |2] shows the time variation of the 
limb width measured on the central detector (CCD#3), for both the direct and reflected 
limb images, during all seven flights. The initial steep portion seen in most flights is during 
ascent, before thermal equilibrium is attained. During Flight 11, in 2009, the direct-image 
limb widened drastically during the float portion of the flight, while the reflected- image limb 
remained steady. The cause of this behavior is unknown; a thermally-induced bending of 
one or more of the BSW surfaces might be to blame, although it would require a strange 
coincidence of surface deformations to adversely affect the direct image and yet leave the 
reflected image intact. While this flight employed a larger balloon, the thermal environment 
for the observations should not have differed substantially from that of the other flights. 
Regardless, because of the FWHM behavior of Flight 11, only observations near the beginning 
and end of the flight will be used in our analysis. Also, in general, it is evident that the limb- 
width profiles of the last four flights are less well-behaved than the corresponding profiles of 
Flights 6, 7, and 8. Possibly the instrument's optical alignment and/or rigidity were altered 
after Flight 8. It should be noted that on-board accelerometers indicate that Flight 8 was 
one of several that ended with particularly hard landings, in excess of 25 Gs. Typically, 
after each flight the detector assembly is removed for photometric calibration. This also 
gives access to the instrument's optical components, which are visually inspected, although 
no intentional adjustments have been made in the optical alignment since the post-Flight 6 
refurbishment at White Sands. While the properties of the instrument are subject to change 
over time, the relative angles of the surfaces in the optically bonded beam splitting wedge 
are expected to remain constant. 



4. Data Analysis 

An SDS flight typically consists of hundreds of thousands of individual exposures, each 
exposure imaging the Sun's limb at ten different position angles, in total, across the seven 
CCDs (seen in Figured]). The 7-/zsec duration exposures are grouped into "cycles" consisting 
of 18 exposures per cycle. The 18 exposures are completed within 0.22 sec and the time 
between successive cycles ranges from roughly 1 to 3 sec, depending on the flight. Pixel data 
from each cycle (7 CCDs times 18 exposures) are stored together in a single file, along with 
housekeeping data for the instrument at the time of observation. The post-flight processing 
pipeline treats a cycle's worth of data at a time, although each exposure within a cycle is, 
for the most part, processed independently. Only in the final step (a correction based on the 
difference in width of opposite limbs of the Sun and its variation throughout the flight) does 
information from one cycle affect the analysis of another. 
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Fig. 2. — Limb edge width as a function of time (relative to apparent solar noon) throughout 
the seven successful SDS flights. The FWHM of the Gaussian fit to the derivative of the 
profile is shown for the direct image (red) and the reflected image (green) on the central CCD. 
The dashed line indicates the theoretical (Rayleigh) resolution limit of the SDS. Observations 
during the initial portion of each flight, prior to thermal equilibrium being reached, typically 
are discarded from the analysis. The sparse data for Flight 10 is the result of an onboard 
recorder failure. The unusual behavior of the direct-image limb width during Flight 11 is 
also noteworthy. 
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The sequence of steps in the processing pipeline is listed below and a detailed description 
of each step is given in the subsections that follow: 

1. adjust photometrically, i.e., flatfield 

2. remove "spike" artifacts 

3. perform preliminary edge detection for all ten limbs 

4. subtract background level and "ghost" image 

5. make final edge determinations for the limbs using the cleaned profiles 

6. correct for bias due to instrumental broadening of the limbs 

7. transform CCD edge coordinates into focal-plane (x, y) 

8. correct (x, y) positions for optical distortion 

9. correct (x, y) positions for atmospheric refraction 

10. fit direct and reflected image measures to two circular arcs 

11. determine the minimum "gap" between direct and reflected limbs 

12. calculate R sun in terms of the BSW angle 

13. correct for Earth/Sun distance 

14. detrend R sun as a function of delta limb width for the entire flight 



4.1. Photometric adjustment 

Prior to most flights, the response of the seven TC101 linear-array CCD detectors 
is calibrated from a series of benchtop observations of a reference sphere. A sequence of 
exposures of varying duration sample the full dynamic range of the detectors. From these 
measures, a quadratic function is fit to the CCD reading as a function of exposure time, on 
a per pixel basis. This is essentially a photometric bias plus nonlinear flatfield correction for 
each detector. Such calibration coefficients were determined prior to flights 6, 7, 8, 9 and 11. 
In some cases, all coefficients were determined in full (Flights 6 and 8) while in other cases 
only the constant (Flights 7 and 8) or constant and linear terms were redetermined (Flight 
11) from the benchtop measurements. 

In the data processing pipeline, each flight's data are reduced using the corresponding 
pre-flight coefficient sets, with the exception of flights 10 and 11; the processing of these two 
flights makes use of the pre-flight 8 coefficients. 

Subsequent to the application of these second-order coefficients, an additional photo- 
metric correction is needed. The TC101 detectors incorporate two amplifiers per detector, 
controlling separately the odd- and even-numbered pixels. The benchtop calibration se- 
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quences cannot account for differences in the bias levels between the two amplifiers as these 
change every time the electronics are powered. Thus, a scalar adjustment is calculated 
based on the empirical differences between neighboring even and odd pixel values, effectively 
putting all pixels on the photometric system of the odd pixels. The seven scalar corrections, 
one per detector, are calculated separately for each exposure within the flight. Figure [3] 
presents a sample profile of the central detector, before (red points) and after (blue points) 
the photometric adjustments are made. What appears to be two curves in the raw data is 
simply the result of the odd/even pixel offset. 

4.2. Spike removal 

Following the flatfield and bias corrections, the observed image profiles contain noise 
in the form of spikes that are one to several pixels wide and of relatively large amplitude. 
These spikes can fool the subsequent edge-detection algorithm, which relies on searching the 
numerical first derivative of the profile. For this reason, the profiles are filtered by examining 
the difference in value between each pixel and the average value of the pixels that fall 4 pixels 
to either side. Any pixel i whose value differs by more than 120 counts has its value replaced 
by the (i + 4, % — 4) pixels' average. Note that the CCD gain is~8 e-/count and the 120-count 
spike threshold is 3 to 4 times the standard deviation of the signal for the portion of the 
CCDs containing the sun's image. In practice, 0.1 % of the pixels are tagged as belonging 
to a spike. The spike-removal procedure effectively eliminates the problem of false edge 
detections. 



4.3. Preliminary edge detection 

We adopt as the conceptual limb edge the inflection point of the profile. This is de- 
termined by finding the maximum (or minimum, depending on whether the edge is a rising 
or falling one) in the first derivative of the profile. In practice, the numerical derivative is 
noisy, creating local maxima and minima that confuse the edge detection routine. Thus, 
in practice, the profile is smoothed before differentiation. We use 20 repeated applications 
of a 3-pt smoothing. This necessary smoothing will unavoidably broaden the limb and, in 
doing so, will shift the inflection point. The shift, or bias, is a consequence of the difference 
in profile shape just before and after the inflection point. This bias will be corrected in 
a subsequent reduction step. Note that even absent the artificial smoothing introduced at 
this step in the reduction, other factors (primarily the instrument response) have already 
broadened the limb significantly. No matter the cause of the broadening, it will result in a 
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Fig. 3. — A sample SDS limb profile, obtained with the central CCD, at various stages of 
the reduction pipeline. The figure shows the raw profile (red); after photometric adjustment 
(blue); after background- level correction (green); and after the "ghost" image is subtracted 
(black). See the text for a description of the different correction steps. 
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bias in the measured inflection point position. The correction is calibrated as a function of 
the overall broadening and thus will correct for the overall bias, i.e., both the unavoidable 
component due to the instrument and that due to the 3-pt smoothing. 

Once smoothed, a simple numerical derivative is obtained, // = (/i+i — /»_i)/2. The 
index of the maximum value of // is identified, and all adjacent points for which /' is greater 
than 0.5 times the maximum value are used to better refine the inflection point location. 
This refinement is by way of a Gaussian fit to the extracted points. The center of the 
Gaussian is taken as the inflection point position. The FWHM of the best-fit Gaussian is 
also stored. See Figure H| for an illustration of the edge-fitting procedure. In this manner, 
preliminary positions for the ten limb edges are found, in terms of pixel position along the 
seven linear-array detectors. 

4.4. Background and ghost image subtraction 

A careful examination of Figure [3] reveals the presence of a low-amplitude "ghost" of the 
main falling limb, lying in the valley between the falling (direct) and rising (reflected) limb 
images and offset from the main image by about 174 pixels. This image is due to a pair of 
internal reflections within a single element of the BSW. (The front and back surfaces of the 
two "flats" from which the BSW is composed are slightly non-parallel, to avoid fringing.) 
It is important to subtract this ghost image from the observed profile since its presence 
can affect the inflection point position, (IPP); primarily that of the falling limb because of 
the direction of the ghost-image offset. A ray-tracing model of the as-designed SDS optics 
predicts an offset of 178 pixels for the ghost image. Empirically, it is found to lie at an offset of 
174 pixels. Its amplitude is a less predictable function of the coatings on the relevant surfaces 
of the BSW; coatings that deteriorate over time. For this reason, the relative amplitude of 
the ghost image is determined for each exposure, while the offset - which depends only on 
the assumed invariant geometry of the wedge - is set at 174 pixels. 

Note that the rising-limb profile also has a contamination from its ghost, although it is 
difficult to detect as the direction of the offset causes it to fall entirely within the high-signal 
portion of that limb. Still, for the sake of consistency, this limb should also have the ghost 
profile removed. 

In order to subtract the ghost-image profile properly, the relative amplitude of the direct 
and reflected limbs must be determined, as well as the relative amplitude of the ghost and 
main images. This is accomplished by measuring the amplitude of the central detector's 
observed profile at key points with respect to the preliminary inflection point positions of 
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Fig. 4. — Illustration of the technique used to determine the inflection point position, IPP. 
The red curve is the smoothed limb profile. Shown in blue is the numerical first derivative of 
the profile. The heavy blue symbols indicate those points used in the least-squares Gaussian 
fit, which is shown by the black curve. The center and FWHM of the fitted Gaussian are 
indicated by the gray line segments. 
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the main limb edges, IPPd and IPPr, as derived in the previous step. 

The observed profiles in all of the SDS detectors are elevated by a background "sky" level 
that is a combination of true scattered light and uncorrected electronic bias, (presumably 
dominated by the latter). The background level is taken to be the minimum value of an 
11-pixel-wide moving average within a conservatively limited range of pixels presumed to 
be free of either (direct or reflected) solar limb and their ghost images. The pixel ranges 
are based on the preliminary edge positions, the preliminary FWHM limb edge width, the 
174-pixel ghost offset, and an additional cushion of 30 pixels. Thus, for detectors 6 and 7, 
which contain only the direct image limb on the low-pixel end of the detector, the signal-free 
pixel range is from the preliminary inflection point position (IPPd) plus the FWHM plus 
the cushion up to 1728 minus the cushion. For detectors 1 and 5, which contain only the 
reflected image limb on the high-pixel end of the detector, the relevant range is from the 
cushion value up to the preliminary IPPr minus the FWHM minus the cushion. Finally, 
for detectors 2, 3 and 4, which contain both direct and reflected image limbs, the signal- free 
pixel range is from the preliminary IPPd plus the 174-pixel ghost offset plus the cushion 
and extending up to the preliminary IPPr minus the FWHM minus the cushion. Figure [3] 
shows a sample profile from detector #3 after subtraction of the background level calculated 
in this manner. 

Having subtracted the background level, it is now possible to calculate the relative 
amplitudes of the various limb images - reflected relative to direct and ghost relative to 
main. This is accomplished by measuring the height of the profile in the central detector 
#3 at four key positions, i.e., just above the "knee" relative to the main direct and main 
reflected limb edges, and at similar pixel offsets relative to the projected positions of their 
respective ghosts. See Figure [5] for a schematic representation of these four key points. A 
simplified geometry for the underlying instrumental limb is shown for the sake of clarity, but 
the strategy is valid for the actual observed limb. 

The offset from the IPP of a limb edge and what is meant by "just above the knee" is 
largely arbitrary; we take its value to be one half the ghost image offset, or A knee =87 pixels. 
Thus, the needed profile amplitudes are 



h = f( x i) where x x = IPPd - A knee 

fa = f( x g) where x g = IPP D + A ghost - A knee 

f2 = f(x 2 ) where x 2 = IPPr + A knee 

fh = f{x h ) where x h = IPPr + A ghost + A knee 
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Fig. 5. — An idealized representation of the direct and reflected solar-limb profiles, as seen by 
the central detector, and their decomposition into "main-" and "ghost-" image contributions. 
The main images are shown in red-hatch shading; the direct-image limb to the left and the 
reflected-image limb to the right. Shaded in gray is the ghost image of the main image, which 
is caused by internal reflection within a single element of the beam splitting wedge. The 
ghost image is offset from the main image by 174 pixels and is greatly reduced in intensity. 
The main and ghost images combine to produce the observed total profile, represented by 
the heavy blue curve. By measuring the amplitude of the observed profile at key points 
relative to the direct- and reflected-limb edges (located at ~560 and ~1060, respectively, 
in the example shown), it is possible to deduce the relative intensity of the direct image to 
that of the reflected image, and the relative intensity of the ghost image to that of the main 
image. 
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Determination of the various fi is made by a linear fit to the observed profile at the 
pixel nearest Xi along with its three neighbors on either side, and then evaluation of the fit 
at %i. With fa, f g , f 2 , and fa in hand, the relative amplitudes are derived in the following 
manner. 

Assume the relative amplitude of the reflected and direct images is T, and the amplitude 
of the ghost profile relative to the main profile is 7. The observed, central-detector profile 
can be thought of as being built up from an intrinsic instrumental solar limb profile that 
has added to it a reflected, offset and scaled copy of itself, and then this combination having 
added to it a ghost copy of itself, with a separate offset and scaling factor. (In reality, there 
are second- and higher-order ghost images also present, but each of these is successively 
diminished by the factor 7 and in practice only the first-order ghost need be considered.) 

For convenience, our formulation actually uses the main reflected image as the base 
template from which the overall profile is to be constructed. Adopting this, there are two 
points in the template profile that are useful to the formulation; call these Ja and fs, as 
shown in Figure The four points measured in the observed profile can be written in terms 
of these two points' amplitudes and the relative scale factors, 

/i = /A/r + 7 / B /r 

f 9 = 7/x/r ^ 

fa = /A 

fh = fB+ if A 

These can be combined to isolate the reflected-to-direct relative scale, T, 

= (f g f g /f2)T 2 + (/1 - f g fa/f 2 )T - fa (2) 
Having found T, the ghost-to-main image relative scale, 7, can be determined, 

7 = (V/ 2 )r. (3) 

Typical values of T and 7 for the SDS flight data are approximately 0.7 and 0.06, 
respectively. These vary somewhat from flight to flight, as shown in Figure [6j presumably, 
this is due to aging of the surface coatings of the beam splitting wedge. 

From the ghost scale factor, 7, and the assumed A ghost = 174 pixel, it is possible to build 
a ghost-free template from the reflected-image portion of CCD#3's profile. Our pipeline 
constructs one template per cycle from the average of the 18 exposures, interpolated and 
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stored at a super-resolution of 0.1 pixels. Knowing the relative scale factors, (T and 7), the 
preliminary main edge positions, (IPPd and JPPr), and adopting a ghost-image offset of 
Aghost = 174 pixels along CCD#3, the template can be used to subtract the ghost image 
from all the detectors' profiles. Note that the ghost-image offset for the outer detectors, 
in pixels, is larger by a factor of 1/ cos(#), where 9 = 27.1° is the angle between the outer 
detectors and the axis of deflection of the beam splitting wedge, which is approximately 
aligned with the inner detectors. Additionally, a pixel scale factor must also be applied to 
the template when used on the outer detectors to account for the slight misalignment of a 
radial for the main image and that of the ghost image. The misalignment is roughly <ft = 1.3°, 
and the correction factor is cos(0). 

Using the above procedure, the ghost image can be subtracted from the observed profile 
in each of the seven detectors. The result for our sample CCD#3 profile is shown in Figure El 

4.5. Final edge determination 

With all pre-processing completed, including the subtraction of the ghost images, the 
inflection point positions (IPPs) are once again determined. This is done as before, using 
the Gaussian fit to the numerical first derivative procedure outlined in Section 14.31 The 
result is ten IPPs per exposure, along with accompanying FWHM measures. 

4.6. Limb-broadening bias 

The shape of the limb edge is asymmetrical about the inflection point. For this reason, 
any mechanism that broadens the limb will impart unequal amounts of influence from either 
side of the IPP and, thus, will shift it. One expects (and finds) a one-to-one correspondence 
between limb width and the net amount of shift or bias. A calibration of this broadening 
bias can be made using a sufficiently accurate synthetic limb profile. The Code for Solar 
Irradiance (COSI, Haberreiter et al. 2008; Shapiro et al. 2010) calculates the solar spectral 
irradiance and has been used to calculate the intrinsic solar limb profile for the passband 
of the SDS. The profile provided by COSI was of sufficient resolution near the limb edge, 
but lacking toward the interior of the solar disk. For this reason, we have mated the COSI 
profile to the analytical solar limb prescription given by Hestroffer & Magnan (1998). The 
resulting synthetic profile is shown in Figure [TJ scaled to nominal SDS pixels. 

Applying the method of Section H~3l (that of Gaussian fitting of the first derivative of the 
profile) provides the IPP and FWHM for the intrinsic limb. The synthetic profile can then 
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be broadened by a specified amount and the fitting procedure repeated, noting the shift in 
IPP from that of the unbroadened profile. The IPP shifts are plotted versus measured limb 
broadening in the inset panel of Figure [7J The shape of the JPP-shift curve is found to be 
insensitive to the functional form of the kernel used to broaden the synthetic profile, being 
characterized uniquely by the FWHM measure of the broadened limb. A best-fit polynomial 
description of this calibration curve is incorporated into the SDS processing pipeline. This 
FWHM-based bias correction is applied to the IPP values of every limb-edge determination. 

4.7. Transformation to Cartesian coordinates 

The IPPs of the ten limb edges are each a one-dimensional measure of position along a 
linear-array CCD. These positions are transformed into focal-plane cartesian coordinates to 
facilitate circle-fitting of the direct and reflected disk images. The CCDs are soldered onto a 
ceramic mounting block, defining their orientations in the focal plane. The endpoints of the 
seven CCDs, i.e., the midlength points of the outer edges of the first and last active pixels, 
were measured in 1988 using the Yale PDS microdensitometer in opaque light-source mode. 
The repeatability of the endpoint measures was found to be ~ 2 microns. (The Yale PDS 
is a sub micron-accuracy machine; difficulty in defining the low-contrast edges of the active 
areas of the SDS detectors was responsible for the 2-micron precision.) 

Note that while 2 microns in the SDS focal plane corresponds to 20 mas, the geometric 
arrangement of the detectors - primarily the fact that the "gap" between direct and reflected 
image edges is measured within individual detectors - lessens the sensitivity of the derived 
solar radius to uncertainties in the detector positions. Monte-Carlo type tests in which 
the detectors were randomly offset and rotated by amounts consistent with the 2-micron 
uncertainty in their endpoint positions yielded just 3.5-mas variation (rms) in the resulting 
radius determinations. This amount, in quadrature, has little influence on what will be 
shown to be an overall 20-mas systematic error estimate for the SDS radius measures. 

Concerned that numerous hard landings after balloon flights might have jarred the 
detector enough to shift the relative positions of the CCDs, the detector assembly was again 
measured in early 2012 using an OGP Avant 600 measuring machine also at Yale. The 
agreement between the 1988 PDS measures and the 2012 OGP ones is at an rms level of 
4 microns, the expected accuracy of the OGP. Thus, it is reasonable to assume that the 
physical locations of the CCDs have remained stable over the course of the SDS balloon 
flights. The 1988 PDS endpoint measures are adopted for our analysis. 

A simple linear interpolation of the IPP, in pixels, and the known (x, y) positions of the 
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CCDs' endpoints provides the necessary transformation. Note that the detector assembly 
is typically unmounted and remounted between flights and an effort is made to align the 
central CCD with the wedge angle plane, such that the direction of offset of the reflected 
image coincides with CCD#3, the central detector. In practice, the alignment is not per- 
fect and the circle-fitting procedure described in Section 14.101 allows us to determine the 
slight misalignment. The actual alignment is taken into consideration when determining the 
separation and minimum gap of the two fitted circles. 

4.8. Optical distortion 

Optical distortion of the as-designed SDS instrument has been determined independently 
using two different ray tracing programs. These agree in their finding that the maximum 
distortion occurs approximately 40 mm from the center of the focal plane, with an amplitude 
of 0.018 mm. As discussed in the description of Step 10, the apparent radius of the direct 
disk image differs from that of the reflected disk image, in general. One possible explanation 
for this asymmetry is an optical field-angle distortion (OFAD) that is not radially symmet- 
ric with respect to the center of the detector plane. We have explored this possibility in 
two separate ways; by adopting the design OFAD but having it offset from the center of 
the detector plane, and by using a ray tracing program to calculate the OFAD under the 
assumption that the Barlow lens is tilted by up to several degrees. 

The first approach allows a reconciliation of the direct and reflected radii for Flights 6 
through 8, given OFAD center offsets of tens of mm. The ratio of radii is more deviant in 
Flights 9 through 12, (R sunD / R sunR =0.992); no amount of center shift can bring the ratio to 
unity. However, the second approach, that of introducing a tilt in the Barlow lens, does allow 
the ratio to be forced to unity for all flights, using three different values of tilt, although the 
tilt angle exceeds 10° for the later flights, an unrealistically large value. Still, it is informative 
that even under such extreme conditions of tilt and OFAD center shifts, the effect on the 
resulting R sun values is slight. 

A series of reductions of the SDS flight data have been made, using a large range of 
OFAD models and corrections. These sets of reductions are used to estimate the possible 
systematic errors in the final R sun values. 
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4.9. Atmospheric refraction 

Even at the SDS float altitude (atmospheric pressure of 3 - 5 mbar), the instrument 
observes through some residual atmosphere and correction must be made for the effect of 
differential refraction over an angle corresponding to the size of the Sun. The formulation 
adopted is based on that of Smart (1979) and is described in detail by Sofia et al. (1994). In 
practice, the correction ranges from about 5 to 25 mas in the vertical direction and remains 
less than 5 mas in the horizontal direction. 



4.10. Circle fitting 

Having been corrected for distortion and refraction, the limb shape is taken to be cir- 
cular; the solar oblateness being ignored at this point. The IPPs from a single exposure 
provide five points along the limb of the direct image and another five points along the limb 
of the reflected image. For each of the two disk images the points spans an arc of about 54°, 
sufficient to define the circles. Each circle, uniquely specified by its radius and the (x, y) 
location of its center, is overdetermined by the five measures. The "geometric" best fit of a 
set of points to a circle minimizes the sum of the squares of the distances of the points from 
the circle. This is a nonlinear least-squares problem, in general, but a numerically adequate 
and computationally fast approximate solution is available to us, given the particular layout 
of the SDS measures. The location of the two outer points, at roughly ±27° relative to the 
central concentration of three CCDs, naturally give the outer points a more prominent role 
in constraining the best-fit circle. We adopt as an approximation to the geometric best-fit 
circle, the mean of the three distinct circles that are defined by the two outer points and 
each of the three inner points, in turn. This procedure assigns a slightly higher weight than 
deserved to the outer points, but, in practice, the random error introduced by this expedient 
approximation is negligible. In this manner, each SDS exposure yields radius and center 
measures for both the direct and reflected disk images. The separation of the two centers 
determines the instantaneous pixel scale for the exposure. The two radii provide a helpful 
diagnostic but do not directly enter into the determination of R sun . Instead, this derives 
from Equation 4 (in Section I4.12p and the measurement of the minimum gap between the 
two disk images. 
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4.11. Minimum gap determination 

The minimum gap between the direct and reflected images is found from two separate 
quadratic functions describing x as a function of y for the three direct-image limb positions 
and the three reflected-image limb positions measured by the central CCDs. The gap is 
first measured in the coordinate system defined by CCD #3 and then corrected for the slight 
angle between this detector and the direction of displacement of the BSW, (as measured by 
the relative displacement of the direct and reflected disk image centers). 

4.12. Calculation of R sun 

Referring to Figure [U the angular half-diameter of the solar disk is 

R sun = W(l - d/D) (4) 

where W is the wedge angle, nominally taken to be 989.47 arcsec, and d and D represent 
the gap and separation. 

4.13. Correction for Earth/Sun distance 

The SDS flights are made during northern- hemisphere fall, at which time the Earth/Sun 
distance is changing rapidly. It is a simple matter to correct the R sun measures within any 
flight to their corresponding values at 1 AU based on the well-known ephemeris of the Earth's 
orbit. The distance at time of observation, in AU, functions as a simple scaling factor for 
the measured R SU n values. 

4.14. Dependence on direct/reflected limb-width difference 

Despite the use of low-thermal expansion materials in the telescope assembly, the re- 
liance on the stability of the BSW reference angle, and the correction of the IPPs for the 
in-flight variation of the measured limb widths, the resultant R sun measures tend to vary 
significantly throughout the cruise portions of the SDS flights. The root cause of this varia- 
tion lies in the differential behavior of the direct and reflected images throughout the flights. 
The SDS operating concept is based on the direct and reflected disk images differing only 
by their being offset by the BSW angle. And yet, the direct and reflected images' observed 
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properties do differ as witnessed by the measured widths of their respective limb images, 
(see Figure |2]). Thus, the final, crucial correction to be made is for the difference in the two 
image systems, as parameterized by the difference in limb widths. Figure [S] shows the strong 
correlation between measured R sun value and AFWHM 3 , the difference in measured FWHM 
of the direct and reflected limbs in the central CCD#3. The measures to be trusted are those 
in which the instrument performs as designed, with the direct and reflected images being 
similar in their properties, i.e., at AFWHM 3 = 0. In practice, instead of limiting ourselves 
to the relatively small number of exposures for which this holds absolutely true, we limit the 
exposures to those with —2" < AFWHM3 < 2" and adjust for the well-determined trends. 

Correcting for the trends seen in Figure [8] with a simple linear or quadratic fit (the 
results are largely insensitive to the form), the final measures of R sun throughout the seven 
flights are shown in Figure The plot scale is identical to that of Figure [TD1 which shows 
the Rsun measures prior to this final correction. 

5. Solar-radius results 

The final SDS measures of the solar half-diameter, after making all necessary corrections 
and adjustments described in the previous section, are listed in Table El Dates of observation 
are expressed in fractional years and the half-diameter measures are given in arcsec. The 
uncertainties listed are la estimates of the combined random and time-dependent systematic 
errorsQ In all cases but one, (that one being Flight 10), the systematic component is dom- 
inant. We estimate the systematic uncertainties from the variation of R sun values resulting 
from limiting cases of assumptions related to the corrections listed in Section 2J primarily 
the treatment of optical distortion. That is, reductions were made in which the assumed 
optical distortion was as designed and with perfect alignment, as well as with distortion 
fields assumed to be offcenter and/or with a Bar low- lens misalignment chosen to force the 
ratio of direct- and reflected-image radii to unity. The standard deviation of the R sun value 
differences was 20 mas and we adopt this as the systematic component of the uncertainties. 
In the case of Flight 10, the small number of recoverable exposures taken at float altitude 
resulted in a formal random uncertainty of up to 20 mas, depending on specifics of the 
data trimming. We have chosen, conservatively, to add the random and systematic error 
contributions for Flight 11 directly, instead of in quadrature, resulting in a total estimated 



1 There is a constant systematic error that we shall ignore, that of the "as-built" reference angle of the 
BSW. We are interested in time variation of the solar diameter and, thus, it is systematic effects that might 
vary with time that are of concern. Whatever its exact value, the wedge angle is assumed constant. 
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uncertainty of 40 mas for this measure. The majority of the data from Flight 11 also had to 
be discarded because of the anomalous behavior of the limb- width during much of the flight, 
but even with these data removed, the formal random error was insignificant relative to the 
estimated 20-mas systematic component. To be specific, the formal random error for Flight 
10 is 20 mas; for Flight 11, it is 0.3 mas; and for the remaining flights it is less than 0.1 mas. 
The systematic component of the total uncertainties listed in Table |3] is assumed to be the 
same for each flight and is estimated to be 20 mas. 

Note that the present results supersede all previous analyses and presentations of SDS 
solar-diameter measures. Those previous studies included only a portion of the components 
of the present analysis; most importantly, they lacked the final adjustment based on the 
difference between direct- and reflected-image limb widths. Another difference of the cur- 
rent analysis is that besides addressing all known optical processes that might affect the 
results, we have explicitly considered the uncertainties corresponding to such processes. As 
a consequence, we judge that these results are robust. 

6. Discussion 

The R sun values from Table |3] are plotted in Figure [HI along with a tracing of the NOAA- 
tabulated sunspot counts over the same period. Given the estimated precision and stability 
of the SDS, a real and significant variation of the photospheric solar radius is demonstrated. 
The alternative possibility that the observed radius variation is due to an unknown solar 
atmospheric effect related to the surface activity is disproven by the fact that the radius 
variations are not in-phase with the activity cycle. 

We note that the SDS is stepped through rotations during each flight, such that the 
diameter measures are made at different solar latitudes. Measurements at all latitudes are 

Table 3: SDS Solar-radius results 



Fit. Epoch R sun @ 1 AU (") 



6 


1992.* 


32 


959.638 ± 0.020 


7 


1994.* 


31 


959.675 ± 0.020 


8 


1995.* 


32 


959.681 ± 0.020 


9 


1996.* 


35 


959.818 ± 0.020 


10 


2001.* 


33 


959.882 ± 0.040 


11 


2009.* 


37 


959.750 ± 0.020 


12 


2011.* 


36 


959.856 ± 0.020 
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included here, thus lessening any possible dependence on low-latitude activity. (It is in the 
determination of oblateness, the subject of a future study, that an effect due to activity is 
likely to be seen, if at all.) By averaging the SDS measures over large fractions of each flight, 
the mean radius over a range of solar latitude is what is actually determined. For most 
of the SDS flights the effective mean solar colatitude measured, $i, is extremely consistent. 
The exceptions to this are Flight 10 (loss of most data due to a recorder failure), Flight 11 
(rejection of anomalous limb profile data), and Flight 12 (restricted rotation sampling due 
to reduced telescope control). The expected flight-to-flight variation in measured radius due 
to this effect, assuming an elliptical limb shape, is 5cos(26i), where 5 is the actual difference 
in equatorial and polar radius. In Flights 6 through 9 the geometric factor, cos(29i), is less 
than 0.1. Assuming a nominal solar oblateness of 8 mas would induce a systematic variation 
from flight to flight of less than 1 mas. Even in Flights 10 through 12, the factor is at most 
0.4, leading to a maximum systematic offset of 3 mas from the uneven sampling in solar 
latitude. This is still insignficant relative to the estimated 20 mas systematic uncertainty of 
each flight's radius determination. We also note that no variation in measured radius as a 
function of SDS telescope rotation angle, or in solar latitude observed, is seen beyond the 
expected ~10 mas variation associated with true oblateness. Thus, we conclude that neither 
solar activity or the solar oblateness can account for the flight-to-flight radius variation 
observed. 

We note that f modes of the Sun can also be used to find a measure of the solar radius. 
Frequencies of solar f modes bear a very simple relationship to stellar structure. They are 
most sensitive to the total mass and radius of the Sun and not to the details of solar structure. 
Solar models constructed with the conventional value of the solar radius usually have f-mode 
frequencies that do not match the Sun. This gave rise to the concept of the "seismic" radius 
of the Sun (Schou et al. 1997, Antia, 1998). The seismic radius basically defines where the 
kinetic energy of the modes is the maximum. The seismic radius of the Sun appears to be 
somewhat smaller than the photospheric radius. More interestingly, helioseismic data show 
that the Sun's seismic radius varies with changes in solar activity (Antia et al. 2000); in the 
sense that the seismic radius decreases with increase in solar activity. By comparison, the 
photospheric radius variations (i.e., those presented here and determined from limb position) 
are much larger and not in phase with the activity cycle. This can be understood by the 
fact that the radius variations are not homologous (Sofia et al. 2005; Lefebvre et al. 2007), 
with different mass shells expanding by different amounts. 

Simultaneous measurements of the photospheric and seismic radii are of great interest 
since they will allow us to calibrate solar models better. Differences and similarities of the 
temporal behavior of these radii should contain information on how and where magnetic 
fields linked to activity change solar structure (Sofia & Li 2001, Li et al. 2003). 
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7. Conclusions 

Accurate, long-term measurement of the diameter of the Sun is an admittedly difficult 
undertaking, frought with issues of instrument stability and calibration. The SDS design 
overcomes these challenges by relying on a beam splitting wedge, constructed using optical 
contact bonding, to provide side-by-side images of the Sun separated by a stable reference 
angle. The dual images transform the relatively large angular measure into a small displace- 
ment in the instrument's focal plane, and also provide a number of diagnostic parameters 
that can be used to monitor and internally calibrate the diameter measurements through- 
out each SDS flight. The SDS data-analysis procedure presented here is complete, allowing 
reliable comparison of results over the decades-long program. 

The balloon-borne SDS experiment has measured the angular size of the Sun seven times 
over the period 1992 to 2011. The half-diameter is found to change over that time by up 
to 200 mas, whereas the estimated uncertainties of the measures, random plus systematic, 
are typically 20 mas. The variation appears as if it might be cyclic, but if so, it is not in 
phase with the solar activity. Thus, the measured variation is not an artifact of observational 
contamination by surface activity. While the SDS measures span 19 years they are sparse, 
making it impossible to say with any certainty that the observed variation is cyclic. The 
results shown in Figure [LTJ call out for additional observations; thus, future (yearly) SDS 
flights are planned. 

The temporal behavior of the Sun's photospheric radius provides a key constraint on 
models of solar structure, particularly with regard to the location, geometry, and evolution 
of subsurface magnetic fields. 

An experiment that has spanned nearly three decades is only possible with the support 
and collaboration of many people and institutions. Beyond the co-authors of this paper, 
we would first like to thank the people who assisted in its conception and design; this 
includes E. Maier, K. Schatten, P. Minott, H-Y. Chiu, and A. Endal. The fabrication of the 
SDS payload could not have been accomplished without D. Silbert. Flight operations were 
helped by many people including D. Pesnell, and W. Hoegy. We are, of course, grateful for 
the substantial financial support provided by a series of grants from NASA and the NSF, 
and more recently, from the G. linger Vetlesen Foundation, and the Brinson Foundation. 
ONES (France) and CNRS (France), which support the PICARD/SODISM mission, also 
provided financial support for the most recent SDS flight and we thank these institutes 
for their contribution. Modeling of the SDS instrument response was greatly assisted by the 
COSI model provided by A. Shapiro and optical modeling by W. van Altena, R. Mendez and 
L. Ramos-Izquierdo. Our work benefited from previous versions of the SDS analysis pipeline 
constructed by J. Zhang, A. Egidi, B. Caccin, and D. Djafer. Finally, we want to acknowledge 
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the outstanding flight support provided by the personnel of the NASA/Columbia Scientific 
Balloon Facility over the years. 
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Fig. 6. — Ratio of the ghost-image amplitude to that of the main image (7) versus the ratio 
of the reflected-image amplitude to that of the direct image (r). The data points are labeled 
by flight number. The "error" bars show the rms of the estimated ratios within each flight. 
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Fig. 7.— Calibration of the limb-broadening bias correction. The black curve displays 
the synthetic limb constructed from the COSI model mated to a Hestroffer & Magnan type 
profile, arbitrarily positioned at an SDS pixel value of 500. The red curves indicate the effect 
of successive smoothing (broadening) of the original profile. Using the method described in 
the text to determine the inflection point position (IPP), i.e., Gaussian fit to the profile's 
first derivative, the shift in IPP is a monotonic function of the limb width, as shown in the 
inset panel. 
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Fig. 8. — Measured solar half-diameter (corrected to 1 AU) as a function of difference in 
width of the direct- and reflected-image limb profiles on the central CCD detector. The 
linear or quadratic fits shown are applied to the flight data as a correction for the observed 
trend. Only data within the fitted range (—2" < AFWHM < 2") are valid for correction; 
data outside this range are discarded. This results in the rejection of a large portion of Flight 
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Fig. 9. — Solar half-diameter (corrected to 1 AU) as a function of time, relative to apparent 
noon, for the seven SDS flights presented here. The flights are color-coded as indicated in 
the legend. With the exception of Flight 10, the data shown represent roughly one in every 
25 SDS exposures throughout each flight. All of the Flight- 10 data that were successfully 
recorded are shown. Compare to Figure [2] for an indication of the portions of each flight 
during which the instrument was equilibrated and stable. 
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Fig. 10. — Same as Figure |9] but after having corrected for the dependence on limb-width 
difference, AFWHM. Only those portions of each flight for which -2" < AFWHM < 2" 
are corrected and retained. This results in much of Flight 11 being discarded, although a 
sufficient fraction of data remains to allow a reliable estimate of iL 1;r) . 
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Fig. 11. — Final values of the SDS-measured solar half-diameter (corrected to 1 AU) as 
a function of time. The error bars indicate la estimates of the combined random and 
systematic uncertainties. For comparison, NOAA-tabulated monthly sunspot numbers for 
the same period are scaled and overplotted. 



